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ABSTRACT 

Inspired by the Kepler planet discoveries, we consider the thermal contraction of planets close to 
their parent star, under the influence of evaporation. The mass-loss rates are based on hydrodynamic 
models of evaporation that include both X-ray and EUV irradiation. We find that only low-mass 
planets with hydrogen envelopes are significantly affected by evaporation, with evaporation being 
able to remove massive hydrogen envelopes inward of ^ 0.1 AU for Neptune-mass objects, while being 
negligible for Jupiter-mass objects. Moreover, most of the evaporation occurs in the first 100 Myr of 
the star's lives when it is more choromospherically active. 

We construct a theoretical population of planets with varying core masses, envelope masses, orbital 
separations, and stellar spectral types, and compare these against the sizes and densities measured 
for low-mass planets, both in the Kepler mission and from radial velocity surveys. This exercise leads 
us to conclude that evaporation is the driving force of evolution for close-in Kepler planets. In fact, 
some 50% of the Kepler planet candidates may have been significantly eroded. Evaporation explains 
two striking correlations observed in these objects: a lack of large radius/low density planets close to 
the stars, and a bimodal distribution in planet sizes with a deficit of planets around 2Rq. Planets 
that have experienced high X-ray exposures are generally smaller than this size, and those with lower 
X-ray exposures are typically larger. The bimodal distribution is uniquely and naturally explained by 
evaporation model, where, depending on their X-ray exposure, close-in planets can either hold on to 
hydrogen envelopes ~ 1% in mass, or be stripped entirely. 

To quantitatively reproduce the observed features, we argue that not only do low-mass Kepler 
planets need to be made of rocky cores overlaid with hydrogen envelopes, but few of them should have 
initial masses above 20M^, and the majority of them should have core masses of a few Earth masses. 
Subject headings: planets and satellites: composition, formation, interiors, physical evolution 



1. INTRODUCTION 

The spectacular success of the Kepler mission 
has yielded thousands of planetary c andidates (e.g. 
IBorucki et ail 120 111 : iBatalha et ail 120131 ). Most of these 
exoplanets are smaller than Neptune (4i?®) and are likely 
to have masses between a few, to tens of earth masses. 
We can now hope to make substantial progress in under- 
standing the origin and evolution of planetary systems by 
studying the interior composition and orbital structure 
of these objects. 

However, well-known degeneracies prevents us from 
resolving their interior composition. A typical Kepler 
planet is likely composed of a dense core and a tenuous 
envelope: the core can be ice-rich, rock-rich or iron-rich; 
the envelope can consist of steam or hydrogen/helium. 
Unlike the case of main-sequence stars, radius measure- 
ment of these systems (all that is possible for most 
Kepler candidates) can not be used to constrain their 
structure. Even in the case where planetary masses are 
measured (via transit-timing- variations, TTV, or radial 
velocity, RV), multiple solutions for the interior struc- 
ture exist. For instance, a Neptune-massed planet, with 
a density of 2 g cm~ 3 , can have either a thin hydro- 
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gen/helium envelope surrounding a dense rocky core, or 
be entirely ice/water d ominated (e.g. [Adams "etalll2008l: 
gers fc Seageril201(3h . 

But this degeneracy can be broken by studying how 
planetary sizes an d densities corre l ate w ith their phys- 
ical environment (jWu fc Lithwickl I2013D . In particu- 
lar, for those planets at close separations from their 
parent star (<0.2 AU), the total received high en- 
ergy irradiation over Gyr time-scales can represent 
a considera ble fraction of their gravitational binding- 
energy fe.g. lLammer et ail 2003: L ccavelier des Etangsl 
120071: iDavis fc Wheatlevi l2009fK If these planets have 
hydrogen-rich envelopes, evaporation of these envelopes 
could markedly reduce their sizes and increase their bulk 
densities, when compared to planets further out. 

Two recent studies have looked for the im- 
prints of ev a porat ion, but reach opposite conclusions. 
iCiardi et all (j2013f ) note that, for pairs of Kepler plan- 
ets that are bigger than Neptune (4i? ffi , these almost 
certainly have hydrogen envelopes), the inner planets 
are more frequently the smaller ones. However, they 
do not observe such a hierarchy in pairs smaller than 
Neptune. Therefore, they concluded that the correla- 
tions may be a result of the planet formation process, 
as opposed to any post-formation pro cess. This con- 
clusion contrasts with that arrived by IWu fc Lithwickl 
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(|2013l) ; in this study, they measure planetary masses in 
a sample of TTV pairs and find that, within pairs, the 
inner planets tend to be denser. Moreover, when con- 
sidering all available mass measurements together, they 
find that planet densities increase with decreasing or- 
bital periods. Therefore, they conclude that the low- 
mass planets observed by Kepler are composed of dense 
rocky cores overlaid with various amounts of hydrogen 
in their envelopes. The same density and radius trends 
are encapsulated in ma ny of the Kepler mult i-planet sys- 
tems, e.g., Kep ler-11 (jLissauer et al.l 12011ft. Kepler -18 
(iCochran et al.ll20TTh. K epler-20 (|Gautier et al.ll2012l ) & 
Kepler-36 (jCarter et al.ll20H) . 

On the theoretical side, there is a large body of liter- 
ature discussing planet evaporatio n, ever since the first 



hot Jupiters were discovered (e.g. ILammer et a! 
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Yellel [200l iTian et all [20051: IMurrav-Clav et al.l 120091: 
Koskinen et all 120101: IQwen feJacksori 120121 ) ~Most of 
these deal with gas-giants, and only a few focus specifi- 
cally on low-mass planets. We summarize two works here 
t hat are of direct relevanc e to our discussion. 

IQwen fc Jack son (2012|) develop realistic models of hy- 
drodynamic evaporation for low-mass planets, includ- 
ing both X-ray and EUV irradiation. They perform 
the calculations on planet models of different densities 
and radii. They identify that only low mass-planets 
(M p < 60 M®) could experience a significant effect from 
evaporation and suggest that evaporation may lead to 
a stability bou ndary below which al l observed planets 
should lie (e.g. IKos kinen et al. 2007). Lacking a ther- 
mal evolution model, they could not make predictions 
on final planet radi us and density. The det ailed evapo- 
ration models allow IQwen fc Jackson] (|2012t ) to calibrate 
the evaporation efficiency and they argue that this pa- 
rameter can vary by orders of magnitude over the stellar 
and planetary lifetimes. Their work also highlight the 
r elative impor t ances of X-ray versus EUV evaporation. 

lLopez et al.l (|2012t) provide the only significant at- 
tempt to study the combined effects of thermal evolution 
and evaporation. Adopting a simplified energy-limited 
formalism for evaporation and a constant efficiency fac- 
tor -the caveats of which are discussed in Section [2J- and 
they also infer the presence of an evaporation threshold. 
They find that since the evaporation time-scale scales 
with planet mass x planet density, that planets with 
mass x density above a critical value wo uld be evaporated 
to bel ow it. But as is demonstrated by IQwen fc Jacksonl 
(|2012fi , a simplistic evaporation model may severely over- 
or under-estimate the mass-loss rates, especially in the 
first stages of the stars' and planets' lives. 

Now with the large Kepler data set, evaporation of 
low-mass planets deserves a better study and a direct 
comparison with the observations. Moreover, it is useful 
to be able to predict the final planet size and density, for 
a given initial model. This allows us to backtrack the 
original planetary structures at formation. To accom- 
plish these goals, we need to improve on previous work, 
in particular, we need to couple realistic calculations of 
evaporation with thermal evolution models. 

We describe our approach in Section 3, after briefly 
reviewing the evaporation theory in Section 2. We then 
compare our theoretical results directly against obser- 
vations in Section 4. Our evaporation theory success- 
fully explains a number of observed facts, establishing 



that evaporation is the driving process of evolution for 
close-in Kepler planets. In fact, some 50% of the cur- 
rently known Kepler planets may have been significantly 
affected by evaporation. Finally, we discuss the caveat 
and implications of our work in Section 5 and conclude 
in Section 6. 

2. PHYSICS OF PLANETARY EVAPORATION 

Planet evaporation can take place through a variety 
of mechanisms: non-thermal escape; thermal Jean's es- 
cape and hydrodynamic escape. Hydrodynamic evapo- 
ration occurs when the density of the heated region is 
sufficiently high so that the gas is collisional even in the 
super-sonic region of the flow. Thus, it is only hydrody- 
namic evaporation that can produce high enough mass- 
loss to affect the structure and sculpt the planet, which 
occurs at high incident fluxes; we focus on this mecha- 
nism. 

As high-energy photons from the star ionize gas in the 
upper envelope (either hydrogen or metals), the newly 
freed electrons heat up the local gas and the atmosphere 
expands. A flow may be initiated that eventually es- 
capes from the gravitational well of the planet. Let the 
efficiency of converting received energy to PdV work be 
77, so the mass-loss rate is simply: 



M = 77 



AGM p a 2 



(1) 



where Lhe is the high energy luminosity (X-rays or 
EUV), Mp is the planet mass, R p is the planet radius 
and a is the separation from the parent star. 

The usual so-called 'energy-limited' approach is 
equivalent of taking r\ t o be an order-unity con 
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stant (e.g. iWatson et all 119811: ILammer et al, 
ILecavelier des Etangsl 120071: lErkaev et al.l 12007^ . Un- 
der such an assumption, the evaporative time-scale is 
Mp/M oc M p p/ Fhe, where Fhe is the high energy flux. 
So Equation [T] lends itse lf to an evaporatio n threshold in 
terms of mass x density (jLopez et al.ll2012l) . 

However, the evaporation efficiency is not always con- 
stant but may vary significantly with planet mass, radius 
and th e ionizing flux (see Section 5). IMurrav-Clav et all 
(|2009f ) demonstrated that, in the case of EUV evapora- 
tion of hot Jupiters, the 'energy limited' approximation 
is only valid at low fluxes. At high fluxes, the mass- 
loss process is controlled by the ionization/recombination 

1 /2 

balance, yielding mass-loss rates that scale as L EUV /a, 
or, the 'efficiency' (77) decreases with flux. This limit (the 
'recombination-dominated' regime) is similar to EUV 
drive n evaporation of gas clumps dBcrtoldi fc McKee 
1 19901) and protoplanetary discs (e.g. Uohnstone et al. 

ll998tlHollenbach et al.ll2000h. 

For X-ray ionization. IQwen fc Jacksonl ([2012D showed 
similarly that the evaporation is not 'energy-limited', but 
rather the 'efficiency' is controlled by line cooling, and is 
a strong function of planet mass. They found that cool- 
ing is most important for high mass Jovian planets, since 
the higher planet escape temperature and larger physi- 
cal size mean flow time-scale is long. Then there is then 
much time to radiate away the received heating when the 
flow is still sub-sonic. For lower mass, Neptune-like plan- 
ets, the escape velocity is much lower and the physical 
sizes smaller, meaning the flow time-scale is shorter, so 
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much less energy is radiated aw ay and the corresponding 
'efficiency' (77) is higher. While lOwen fc Jackson] (|20f 2fi 
found that the X-ray evaporation gives a similar flux- 
scaling as given in Equation [I] the flow is not close to 
being 'energy-limited' (where PdV work is the dominant 
energy loss channel), and only begins to approach an 
energy-limited case at low masses < 3 M®. Furthermore, 
Owen & Jackson noted that unlike the EUV case, since 
the X-ray driven mass-loss scaled as Lx/a 2 at high and 
low fluxes there is no transition from line-cooling limited 
to energy limited at low X-ray fluxes. 

The ionizing flux from a main-sequence star varies by 
orders of magnitude through its life. In addition, the 
size of a planet also changes over time, under both ther- 
mal contraction and mass-loss, along with the mass of 
the planet decreasing. With r\ being a function of planet 
mass and radius, as well as the ionizing flux (see Fig- 
ure [TT|), one should not adopt a constant 77 (as in Lopez 
et al. 2012) in tracing out the evaporation history of 
the planet. It is important that the correct prescription 
of evaporation is used to follow the planetary evolution, 
particularly for low mass planets where evaporation can 
significantly affect the evolution. This is essential if one 
wants to make inferences about the initial state of these 
planets. 

The last issue worth attention is the nature of the ion- 
izing flux. Which source of high energy luminosity is 
driving the evaporation has only been tackled recently. 
Som e authors have purely used E UV flux in their studies 
(e.g. lLecavelier des Etan gs 2007) , wh ile others have used 
only t he X-rays ([Jackson et al.N2012D . [O wen fc Jackson! 



(|2012f) solved the flow problem including both X-ray and 
EUV irradiation. They found that the position of the 
transition from sub-sonic to super-sonic flow (the 'sonic 
point'), relative to the respective ionization fronts, deter- 
mines whether X-rays or EUV is driving the mass-loss. 
At a similar ionizing energy flux, they found that X-ray 
drives the mass-loss when the flux is high, and EUV dom- 
inates the mass-loss when the flux is low. For the first 
~ 100 Myrs of their lives, main-sequence stars have X-ray 
luminosities that reach up to ~ 10~ 3 of their bolomct- 
ric luminosities; this fra ction falls off with time roughly 
as 1/t as the stars age (jGiidell 12004 iRibas et all 120051: 
Uackson et aTll201 2[ ). 

3. EVOLUTIONARY MODELS OF PLANETARY 
EVAPORATION 

In order to follow the evolution of a close-in planet, 
we need to include the effects of thermal cooling of the 
planet, irradiation of its upper atmosphere and mass- loss 
in the form of hy drodynamic evaporation. Previously, 
lLopez et all (120121) have address ed thi s problem by cou- 
pling the iFortnev et al.l (|2007bl ) and iNettelmann et al.l 
(|2011l) planetary structure models to a simplistic energy- 
limited estimate for the mass-loss. We aim to achieve a 
similar goal here, but we will base our mass-loss rate on 
a realistic calcul ation that include both the X-ray and 
EUV radiation (jOwen &: Jackson] 120 12T). We make use 
of th e mesa stellar evolution code (jPaxton et aD I2011L 
2013) to model the thermal evolution and couple it to 
the Owen & Jackson mass-loss calculations. 

3.1. Method 



The general purpose mesa code provides the frame- 
work to simulate planetary structure and evolution. Our 
planet models have solid inner cores, experience irradia- 
tion (under the two stream approximation, e.g. iGuillotl 
I2010T ) . and evaporation. 

To include evaporation, we tabulate the mass-loss rates 
as functions of ionizing fl ux, planet mass a n d pla net ra- 
dius, using the results in lOwen fc Jackson] (|2012[ ). The 
tabulated results span a range in planet mass from 1M® 
to 3Mj, and in planet radius from li? ffi to ~ 0.01 AU 
(roughly the Hill radius at 0.1 AU for our most mas- 
sive planet). For t he ionizing X-ray flu x, we adopt the 
observed values in (Uackson et al.ll2012f ) , which are func- 
tions of both stellar spectral type and age. In general, 
the X-ray luminosity saturates at <~ 10~ 3 of the stellar 
bolometric luminosity during the first ~ 100 Myrs of the 
star's life, and decays approximately as 1/t af terwards. 
Furthermore, following lOwen fc Jackson] (|2012l ). we take 
the EUV luminosity to follow the same time evolution of 
the X-rays. 

Since convective transport is the dominant heat trans- 
port mechanism in the planet's envelope, the thermal 
structure of the envelope is set to be initially adiabatic, 
in the absence of irradiation. The radius of the planet is 
defined to be at a radius where the optical depth to the 
incoming bolometric radiation is 2/3, and is typically 
around miliBar (for younger and lower mass planets) 
to Bar (for older and more massive planets) pressures 
for the planets considered here. This planet radius is 
also taken as the input radius in the evaporation model. 
Since the atmosphere's underlying scale height is typi- 
cally small compared to the planet radius, such an ap- 
proximation is accurate. But it may break down for the 
closest in planets at the earliest times, where there may 
be a ~ 10% difference in a planet's optical photosphere 
and the base of the evaporative flow. However, given 
that planets quickly contract through such a phase, the 
effect will be negligible when integrated over the Gyrs of 
planetary evolution. 

Our planets are composed of two separate regions: an 
envelope of Hydrogen/Helium with a mctalicity at solar 
abundances (as used for the evaporation models), and a 
solid core. Given the large range of possible core compo- 
sitions, as a starting point we focus on a pure rock core, 
using t he mass-radius profile provided by IFortnev et al.l 
(|2007al|bT ). with core masses (M c ) of 6.5, 7.5, 10, 12.5 
& 15 M® adopted. As current planet formation mod- 
els are unable to provide a good handle on the initial 
thermal properties of formed planets (e.g. 'hot' or 'cold' 
start) we consider initial models with a wide range of 
initial radii. This initial radii (or more correctly en- 
tropy) is parametrised in terms of an initial cooling time 
(£cooi), which we define as the ratio of a planet's ini- 
tial internal energy to a planet's initial luminosity. We 
could of course parametrise the initial entropy in terms 
of some other variables; however, the initial cooling time 
is perhaps the most interesting as this be compared to 



protoplanetary discs lif etime s ~3 Myr (e.g. lHaisch et al 
' 200lMHernandez et al.ll2007tlOwen et alJl201lHArmitag" 
20111 ). We choose values of i cool (computed in the ab- 



sence of irradiation) in the range 3 x 10 6 — 10 8 years to 
span a range of 'hot' start and 'cold' start model, similar 
to those values chosen by Lopez et al. (2012). We note 
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Fig. 1. — The top panel shows the radius evolution of Jupiter like 
planets as a function of time since gas disc clearing, at a very close 
separation to their parent star (~ 0.025 AU); the bottom panel 
shows the mass evolution of this planet. The solid line represents a 
planet with an initially high entropy with an initial cooling time of 
10 6 years, the dashed line has an initial cooling time of 10 7 years, 
and the dotted line has an initially low entropy with an initial 
cooling time of 10 s years. 

that with cooling times < 3 x 10 6 years most low mass 
planets at separations < 0.1 AU have initial radii larger 
than their Hill radii, and they cannot be considered hy- 
drostatically bound objects. 

For each cooling time and each core mass, we construct 
40 models with a range of envelope masses from 3 Mj 
to a few percent of the core mass, logarithmicly spaced 
in envelope mass. The planets are then evolved forward 
in time, under the influence of evaporation and irradi- 
ation for 10 Gyrs or until the entire envelope is lostQ. 
Our integrations begin at 3 Myrs, a time at which the 
protoplanetary disc clears and the planet is fully exposed 
to the X-ray and EUV irradiation. In general the plan- 
ets evaporation begins in the X-ray driven phase, and 
may switch to EUV driven at some late time. Where the 
planets' final mass is set by a few 100 Myrs. 

3.2. Jupiter-like planets 

Jupiter mass planets close to their central stars repre- 
sent the case where there is direct ob servational evidence 
(e.g. IVidal-Madiar et al.l 12001 120041 ) of evaporation oc- 
curring. Thus, it is worth investigating whether their 
are any evolutionary consequences for the evaporation 
of high-mass planets. For example, iBaraffe et al.l (|2004f ) 
noted that if the evaporation rate was high enough, so 
that the evaporation time (~ M p /M) became compara- 
ble to the thermal time of the planet's envelope, then 
a Jupiter mass planet could loose its entire envelope 

Thus, we ignore any subsequent mass-loss from the core due 
to sublimation, which may h appen at the highest equilibriu m tem- 
peratures T eq > 2000 K (e.g. IPerez-Becker fc Chiandl2Tjl3h . 



rapidly. Such an inference was based upon the rather 
unrealistic assumption of 100% mass-loss efficiency. In 
reality, at high masses the mass-loss rates never reach 
such high values (see Owen & Jackson 2012 for a more 
detailed discussion) . To illustrate this we show the evolu- 
tion of a Jupiter-like planet at a very close separation of 
0.025 AU around a solar type star in Figure [1] We plot 
the radius and mass evolution of planets with a 15M0 
rock core and cooling times ranging from the very short 
(10 6 years - solid) to the very long (10 8 years- dotted). 
Figured] clearly shows that evaporation is unable to effect 
the planets' evolution and planets with very different ini- 
tial cooling times end up on almost identical evolutionary 
tracts at late times. With the mass-change in the planets 
at the < 1% level. 

Therefore, while Jupiter mass planets currently pro- 
vided the best opportunity for actually studying the hy- 
drodynamics of evaporation by directly probing the flow, 
they do not provide a good laboratory for studying the 
evolutionary consequences of evaporation and we must 
go towards lower mass planets where the evolutionary 
effect will be more pronounced. 

3.3. Low-mass planets 

Unlike the Jupiter-type planets discussed above, sev- 
eral authors have suggested that the effects of evap- 
oration will be more prominent around lower mass 
planets (e.g. IBaraffe et~aTl 120081 I Jackson et all 120121 : 
lOwen fc Jackson! I2012t lLopez et al.1 l2012h . At lower 
masses evaporation can begin to sculpt the planet pop- 
ulation, removing significant amounts of a planet's en- 
velope during its lifetime. To investigate this we show 
the evolution of a low mass planet in Figure [2j where we 
consider the evolution of a 'standard' model which is an 
initially 20 M ffi planet with a 12.5 Mq core. Figure [2] 
shows the evolution of the 'standard' model for the full 
range of initial cooling times considered. The planets are 
at a close separation of ^0.05 AU and have equilibrium 
temperatures of 1300K. Where we define the equilibrium 
temperature as the black body temperature at a given 
separation, i.e.: 




where T* and i?* are the parent star's temperature and 
radius respectfully 

The panels in Figure [2] show the qualitative features 
of the evolution of low-mass planets. In particular, plan- 
ets with higher initial entropies have initially larger radii 
and therefore loose more mass. The result is that plan- 
ets with shorter initial cooling times end up with, smaller 
radii and higher densities than planets with longer ini- 
tial cooling times. At low masses evaporation plays a 
strong role in planetary evolution, which in the case of 
our 'standard' model makes the planets a factor of ~ 2 
smaller and a factor of ~ 4 denser compared to the same 
planet that is not undergoing evaporation. In this case 
an initial envelope containing ~ 40% of the original mass 
is evaporated down to one containing only few percent 
of the total mass, and in the most extreme case (for the 
planet with an initial cooling time of 3 x 10 6 years) evap- 
oration leaves a planet with 0.5% of the total mass in a 
Hydrogen/Helium envelope. 
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Fig. 2. — The evolution of a 20 Mj planet with a 12.5 Mq core at a separation of 0.05 AU is plotted, with four initial cooling times: 
3 X 10 6 years (dashed); 1 X 10 7 years (solid); 3 X 10 7 years (dotted); 1 X 10 s years (dot-dashed). In the upper left we show the radius 
evolution, the density evolution is shown in the upper right, the mass evolution in the bottom left and in the bottom right we show the 
time evolution of M X age which indicates the times at which most mass-loss is occurring, with the thin continuation to these lines show 
the point at which the evaporation switches from X-ray driven to EUV driven. The thin, straight, solid line show the time at which the 
X-rays are saturated until, and the final properties of the same planet not undergoing evaporation are shown by the small arrows. 



The bottom right hand panel of Figure[2]shows the evo- 
lution of M x age, which indicates when the mass-loss is 
most significant. This plot shows that in all cases the 
mass-loss is most important at roughly the point where 
the X-ray luminosity begins to decline. This is easy to 
understand as at early times the planets are large and 
the fluxes are high, so the planet can absorb a significant 
amount of high energy radiation. Once the X-rays begin 
to decline the planet evolution is affected less by evapo- 
ration and once the evaporation switches to EUV driven 
the final planet properties have already been 'frozen' in, 
similar to the results from the previous models calculated 
by Owen & Jackson (2012) and Lopez et al. (2012). This 
is a rather generic feature of all the evolutionary models, 
where the saturation time-scale for the X-rays sets the 
length of time which evaporation is important in driving 
the planets' evolution, and that there is very little change 
in the planets' mass from 100 Myrs to 10 Gyrs. 

4. POPULATION STUDY AND COMPARISON AGAINST 
OBSERVATIONS 

We have shown that, for low mass planets close to their 
parent star, a significant amount of the planet's gaseous 
envelope, or even an entire atmosphere, can be removed 
in Gyr time-scales. Given the prominent role evaporation 
can play, we can use this evolution to make inferences 
about the initial structure of observed close-in planets, 
and provide clues as to their formation. In the following, 
we provide comparisons between our theoretical mod- 
els and observations, mostly using data from the Kepler 
mission. The effects of evaporation are clearly visible in 
data. In fact, the evaporation model naturally explains 



a number of remarkable correlations seen in the Kepler 
data. 

In Figure [31 we plot the final planet mass and radius as 
a function of separation for planets with initial proper- 
ties which vary from our 'standard' model, where central 
star has 1 M in all cases. At a given planet mass, a 
lower core mass results in initially larger planets, which 
means that the planet can absorb a larger fraction of the 
X-ray flux, driving stronger evaporative flows. So plan- 
ets with smaller cores may lose the entire envelope at a 
larger separations than those with larger cores (see Fig- 
ure [7]) . The initial cooling time shows the same effect 
as discussed above, with shorter cooling times resulting 
in more mass-loss due to the initially larger radii; how- 
ever, this effect is much less important than the effect 
of initial core mass. Finally, evaporation can drive con- 
vergent evolution. For instance, planets at < Q.05AU 
that had initial 15 M ffi cores, but envelopes of 2 & 5 
M ffi end with somewhat identical structures This means 
that, at least in some cases, it may be difficult to retrieve 
the planets' initial structure, particularly at low envelope 
mass-fractions. 

4.1. Upper Envelope in Planet Radius 

The Kepler transiting-planet catalogue is the most ex- 
tensive catalogue of planets on close orbits to their parent 
star; for better statistics, we use the Kepler object of in- 
terest (KOI) catalogue, which lists the radii (not mass) 
of planet candidates. Most of the KOIs have not been 
confirmed as planets, and there is a certain percentage of 
false positives. To minimize contamination, we choose to 
consider only KOIs that have been identified as being in 
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Fig. 3. — Final planet radius (top panel) and mass (top panel) as 
a function of separation, for models with different initial core mass 
(labelled as M c ), total mass (labelled as M p ) and initial cooling 
time (thin lines for 3 X 10 6 yrs and thick lines for 10 8 yrs, longer 
cooling time indicates lower internal entropy). The central star is 
sun-like. The straight dotted lines show the masses and radii of 
the rocky cores alone. Planets inward of ~ 0.1 AU experience sig- 
nificant mass-loss and can retain only a fraction of their primordial 
atmospheres (depending on the core mass), with the hottest ones 
exposed to bare cores. As a result, the maximum radii of these 
planets decrease inward. 

multiple systems, where the signifi cance of planetary na - 
ture is considerably higher (> 95% iLissauer et aL1l2012[ ). 

Plotting planet radius against separation for planets 
in multi-planet systems around solar-type stars (T* = 
5200 — 6200 K) we notice an upper-envelope in planet 
siz e that rises with sep arati on. This has also been n oticed 
bv lCiardi et all (f2013h and IWu fc Lithwickl (pOll . The 
lack of large planets at small separations is statistically 
significant: we cut the sample into 4 regions: large and 
small planets with a division at R p — 5 and hot 
and cold planets with a division at T cq = 1000K, then 
comparing the ratio of large cold planets to the ratio of 
small cold planets. Where one would expect to observe ~ 
20 hot large planets if these are the same populations, we 
find only 3. This upper envelope is naturally explained 
by evaporation: low density planets can not survive in 
the environment of high ionizing flux. Quantitatively, 
evaporation of 20Af ffi planets with rocky cores of masses 
10-15 Mq provides a good fit to the upper envelope, as 
is shown in Figure H] 

There is another implication to this agreement. If 
there were a significant population of planets with ini- 
tial masses higher than, say, 30M e (dotted line in Fig- 
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Fig. 4. — The upper envelope of planet sizes as a function of 
equilibrium temperature . The open circles show the Kepler Object 
of Interests that are around solar type stars (T* = 5200 — 6200 K) 
and are in multiple transiting systems. The black curves are the 
theoretical final radii for planet models with initial mass 20 
and 10 (dashed), 12.5 (solid) & 15 M® (dot-dashed) of rocky cores. 
The dotted line are for ~ 30 planets with 12.5 Mq cores. All 
models here have an initial cooling time of 10 7 years. 

ure[4]), we would not expect to observe the same upper 
envelope. These more massive planets can hold on to 
their atmospheres, much like the hot Jupiters can (see 
Section 13. 2[) , and would populate the upper left region 
in Figure |4j The absence of these more massive plan- 
ets is interesting: perhaps not coincidentally, as 20 — 30 
M0 roughly corresponds to gap-opening mass in this re- 
gion, suggesting that there is ceiling to how much gas 
the planets can accrete in this neighbourhood. The core 
mass also somewhat affects the upper envelope. We find 
that models with roughly half of the total mass in the 
rocky cores b est reproduce the data, consistent with the 
conclusion in IWu fc Lithwic k (2013). In contrast, the 
planet's initial cooling time makes little difference to the 
results. 

4.1.1. Dependence on Stellar Spectral Types 

The mass of the parent star is an important consider- 
ation. The parent stars direct influences on the evapo- 
ration, through its gravity is small. However, the total 
received X-ray flux at a fixed bolometric flux (so fixed 
equilibrium temperature) can vary greatly, with late-type 
stars being significantly more X-ray luminous when in- 
tegrated over Gyr time-sc ales (|Giidell l200l iGiidel et all 
120071 iJackson et alJf20r2T ) . We demonstrate this by con- 
sidering the evolution of the 'standard' planet discussed 
above (M p = 20 M e , M c = 12.5 M e ), at a fixed equi- 
librium temperature (1300K) around a later-type (0.5 
M Q ) and earlier-type (1.5 M Q ) star compared to a planet 
around a solar type star, all with initial cooling times of 
10 7 years. 

The evolution of these planets, is shown in Figure [5j 
where the radius evolution is in the top panel and the 
mass evolution in the bottom panel. 

Naively, one would expect a similar evolution as the 
bolometric flux revived is identical in all cases. However, 
the variation of X-ray luminosity with stellar mass results 
in qualitatively different evolutionary paths for the plan- 
ets, with order-unity differences in both the final planet 
mass and final planet radius. The planet around the 
0.5 M star has had its envelope completely removed, 
whereas the planet around the 1.5 M Q planet still has a 
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Fig. 5. — The panels show the radius (top panel) and mass (bot- 
tom panel) evolution of a planet (20 Mq) with an equilibrium 
temperature of 1300K (~0.05 AU around a solar type star) and 
12.5 Mq rock core, around a 1.5 Mq star (dot-dashed), 1 Mq star 
(solid) & 0.5Mq star (dashed). The solid lines are identical to the 
solid lines in Figure [2] 

~ 3 envelope remaining after 10 Gyrs. 

We can go further and compare our evaporation thresh- 
old found above for solar type stars to KOIs around other 
types of host stars. Lower mass stars (e.g., M-dwarfs) 
have higher X-ray flux, when compared to their bolomct- 
ric luminosities. They also remain choromospherically 
active for longer periods. Figure [5] shows that, indeed, at 
the same equilibrium temperature (measuring the bolo- 
metric flux), planets around lower mass stars are more 
severely evaporated. 

For a quantitative comparison, we determine the the- 
oretical evaporation threshold by extracting a linear re- 
lation between the maximum radii and equilibrium tem- 
perature, for planets with an initial mass of 20 and 
a core mass in the range 10 — 15 M®, that are orbiting 
around a solar-type star and have equilibrium tempera- 
ture in the range of 500 — 2000 K. Since the total mass 
loss roughly scales linearly with the integrated X-ray flux, 
we expect the same radius threshold to apply to planets 
around all spectral types, when we arrange them by their 
X-ray exposure. This is shown in Figure [6l where the 
planets are separated by the spectral type of their parent 
star0 In contrast, we show the same planet radii plotted 
against their bolometric exposure. Planets satisfy the 
same evaporation threshold only when one considers X- 
ray exposure. This argues that the ionizing flux, not the 
bolometric flux, is what determines the upper envelope. 

4.2. Distribution of Radius 

2 Since the evolution of the X-ray luminosity is poorly known 
for stars < 0.45 Mq , we approximate the X-ray evolution of these 
stars with that of a 0.45 Mq star. 
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Fig. 6. — The observed radii of KOIs are plotted against their 
bolometric exposure (bolometric flux received at surface integrated 
over 10 Gyrs, le ft panels), and X-ray exposure (same but for X-ray 
exposure, using rjackson et al.1120121 values, right panels). The ob- 
jects are roughly separated into those around A/F stars (top pan- 
els), G/K stars (middle- panels) and M-dwarfs (lower panels, the 
red squares stand for the system KOI 961). KOI multis are shown 
as open circles while KOI singles as small dots. As predicted by 
our models of x-ray evaporation, planets around stars of different 
spectral types, while having very different bolometric exposure, 
have roughly the same maximum sizes at a given x-ray exposure. 
The solid lines show the theoretical evaporation threshold derived 
from a linear fit to the black evaporation curves shown in Figure [4] 
The right panels demonstrate that planets around different types 
of stars satisfy the same maximum radius/X-ray exposure relation 
as that found around solar-type stars. One cannot perform the 
same exercise in the left panels. 
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Fig. 7. — The final radii of planets with an initial mass < 20Mq, 
as a function of separation from a sun-like star. Different colours 
represent different core masses with blue standing for 15 Mq core, 
red for 12.5 Mq core, black 10 Mq, magenta 7.5 Mq & green 
6.5 Mq. We consider atmospheres as low in mass as one percent 
of the core mass. For all models, there is a critical separation 
within which the planets are evaporated down to bare cores. This 
separation is larger for models with smaller core masses. 
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Our analysis in the previous sections suggest that the 
observed radius cut-off in KOIs is related to the fact that 
the most massive KOIs have masses not much exceeding 
2OM and that their core masses are roughly half of their 
total masses. Here, we investigate the nature of all KOIs 
by studying the overall radius distribution. 

In Figure we present the final radii of multiple se- 
quences of planet models with initial cooling times in the 
range 3 x 10 6 — 10 8 years. These have core masses from 
6.5 M ffl to 15 M®, and atmosphere masses from approx- 
imately one percent of the core mass to much larger val- 
ues. The total planet masses are restricted to < 20Af®. 
We do not consider atmospheres with masses below one 
percent, motivated by the discussion below. The radii 
are evaluated after 10 Gyrs of orbiting around a sun-like 
star, though th e values differ little if we instead eyaluat e 
at 1 Gyrs (e.g. lLopez et al.ll2012t lLammer et aT1l2013h . 
Figure [9] shows the corresponding planet densities as a 
function of separation and is discussed in Section 14.31 

The overall population show the general feature noted 
previously, that the radius decreases with decreasing sep- 
aration. One particularly interesting feature that ap- 
pears is a gap in radius between planets that have gaseous 
envelopes, and those that are bare cores. For instance, 
for the 6.5 M® core models, inward of <~ 0.1 AU, all 
planets have their atmospheres stripped away with their 
final sizes reflecting that of their naked rocky cores; while 
outside ^0.1 AU, planets can retain atmospheres that 
are at least a fraction of a percent in mass, consequently 
they have sizes that arc markedly larger. This bifurca- 
tion generates a gap in planet radius. The orbital sep- 
aration at which this gap appears is smaller for planets 
that have bigger cores and stronger surface gravity. How- 
ever, inside ~ 0.03AU there are no surviving planets with 
gaseous envelopes. 

The origin for this gap is easy to understand. As plan- 
ets lose their hydrogen atmosphere, they become increas- 
ingly compact and dense, which reduces the mass-loss 
rate. However, there is also less hydrogen to lose. So 
for planets inside the critical separation, the loss is to- 
tal. While for planets just outside the critical separation, 
there is a minimum atmosphere mass the planets need to 
avoid complete stripping. Any thinner atmosphere will 
be easily eroded. This mass is roughly 1% of the total 
mass and corresponds to roughly an order unity modifi- 
cation to the planet radius. So we expect to see a gap in 
planet size. The presence of such a gap naturally arises 
in evaporation models, but is not expected in other sce- 
narios, for example planet mergers. 

In Figure [HI we demonstrate that the radius distri- 
bution of KOI multi-planet systems is bi-modal. Most 
planets have sizes ~ 1.5-R® or 2.5i?®, with a deficit of 
objects at radius ~ 2i?®, indicating the presence of a 
gap. To confirm this suspicion, we divide the objects 
by their X-ray exposure into a high X-ray group (corre- 
sponding to < 0.1 AU around a solar- type star) and a 
low X-ray group. Strikingly, objects with a high X-ray 
exposure mostly have sizes below the gap, at ~ 1.5-R®, 
while objects with a low X-ray exposure typically have 
sizes above the gap, at ~ 2.5i?©. This argues that the 
deficit at 2i?® is physical and is associated with X-ray 
exposure. 

We believe such a connection is real and is not related 
to various selection bias in the KOI catalogue. The same 
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Fig. 8. — The correlation between radii of KOI and their X- 
ray exposures. In the left panel, we show the radii of KOIs in 
multi-planet systems (in open circles) versus their integrated X-ray 
exposure. The division between the red and the green population 
corresponds to a distance of 0.1 AU around a sun-like star. The 
black histogram in the right panel is the size distribution of all KOI 
multi-planet systems, while green that of the low X-ray group and 
red the high X-ray group. The red and green objects have distinctly 
different peak sizes, as is expected from our X-ray evaporation 
models (Figure 0. While the green objects (with a peak size at 
~ 2.5-Re) manage to hold on to some of the hydrogen envelopes 
they were born with, most of the red objects have been stripped 
down to naked cores (with a peak size at ~ 1.5-R®). The bimodal 
size distribution of KOls is naturally explained by their mass-loss 
history (see text). 

bi-modal behaviour appears when we consider only KOI 
singles, the full KOI catalogue, or when we include only 
bright KOI targets, or only dim KOI targets. We further 
note that although the values for stellar radii (and there- 
fore inferrred planet radii) may be erroneous for some 
candidates, such errors are expected to diffuse any fea- 
ture rather than to generate it. Better stellar determina- 
tions ma y further sharpen the b imodal feature (see recent 
work by [Morton fc Swift] [20131 ) 

The presence of this gap provides strong evidence that 
evaporation not only sculpts the upper envelope of planet 
sizes, it is also driving the evolution of the majority of 
KOI objects, with ~ 50% of Kepler planet candidates 
having experienced significant mass-loss during their life- 
times. Comparing the gap location (0.1 AU around sun- 
like star) against our theoretical calculations, we suggest 
that the planet population have predominately low mass 
cores ~ 6 M® , and that most started out their lives with 
Hydrogen/Helium envelopes of at least a few percent in 
mass. 

Figure [7] illustrates that the gap radius, as well as the 
separation at which this gap appears, are direct probes 
of the core properties. An improved investigation on core 
composition and mass should be conducted when planet 
radii are better determined. 

4.3. Comparison of Planet Density 

While the KOI catalogue only allows comparison of 
planet radius, the small but growing sample of low-mass 
planets with measured masses also provides another im- 
portant comparison: planet density. In Figure [5J we il- 
lustrate planet densities resulting from our integration. 
In the density-separation plane, the upper envelope in 
planet size is now translated into a lower envelope in 
planet density. There is a gap in planet density, similar 
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Fig. 9. — Same as Figure [7] but showing final planet density. 
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Fig. 10. — Observed planet density plotted against orbital separa- 
tion, scaled to have the same X-ray exposure as around a sun-type 
star. The planet sample includes both th ose from radial ve locity 
surveys (red filled circles, as compiled by [Weiss ct al. 2013|) and 
Kepler pl anets with masses confirmed by the Kepler team (com- 
pilcd from Lissaucr ct al. 2011; Cochran ct al. 2011; Gauticr ct al. 
2012; Carter ct al. 2012) or calculated from transit timing varia- 
tions (TTV) bv lWu LithwTckl j20T3l) . There are large error bars 
to many of the TTV masses. The shaded region is the disallowed 
region as shown in Figurc[9] and the dashed lines show the densities 
of bare rocky cores for a range of masses. 

to that in radius. However, if there is a planet popula- 
tion with low mass cores (1 — 3Mgg), these low density 
cores may partially fill in the gap. 

To compare, in Figure [TU] we plot the measured densi- 
ties of a list of low masses planets against their scaled sep- 
arations. We scale the actual planet separation by their 
respective X-ray exposure, though since most host stars 
are solar-type, this correction is typically small (< 10%). 
The theoretically disallowed low-density region is shown. 
Most of the observed densities avoid this forbidden re- 
gion, providing strong evidence for sculpting by evapo- 
ration. 

Current density measurements (especially those using 
transit-timing fittings) contain large uncertainties. This 
prevents us from making more quantitative comparisons 
at the moment. In particular, we could not discern the 
density gap as predicted by model ca lculations. 

The particular case of Kepler-36bc (|Carter et al.ll2012T) 



is worth commenting. The two planets have nearly iden- 
tical separations (0.115 AU and 0.128 AU) but largely 
discrepant densities, 6.8 & 0.86 g cm -3 , respectively. At 
their present orbits, the minimum core masses for the 
two planets to retain their envelopes is ~ 6.5 M ffi (see 
Figure [9j. The measured masses of the two planets arc 
4.45 M ffi and 8.08 M ffi , respectively, naturally explain- 
ing the diversity in their structure. More systems like 
Kepler-36 will be able to provide strong constraints on 
the nature and strength of evaporation in close-in plan- 
etary systems. 

5. DISCUSSIONS 

We discuss some of the caveats and limitations of 
the presented calculations and how they bear on our 
inferences about the observed planet population. The 
two most important assumptions concern the use of 
a two-layer planetary model (a rock core plus a Hy- 
drogen/Helium envelop e), and the adoption of the 
IQwen fc Jacksonl (|2012f ) evaporation model. 

5.1. Variations in the planetary structure 

Low-mass planets could contain volatile-rich atmo- 
spheres, they could also contain a significant amount 
of iron and/or ice in their cores, ()Adams et al.l 120081 : 
iRogers fc Seagerll2010D . Our simple two- layer model of 
rocky cores plus hydrogen envelope has been successful 
in explaining the observations. But what about these 
other possibilities? Can we exclude them based on cur- 
rent observations? 

We can exclude the possibility that the dominant 
primordial atmospheres of these low-mass planets are 
very rich in volatiles. Our arguments b elow run sim- 
ilarly to those given in IWu fc Lithwicl (120131) but arc 
more informed by our detailed modelling and by the 
physics of evaporation. If one considers a water-rich 
envelope, for instance, evaporation will not proceed as 
is described here. First, water molecules have to be 
photo-dissociated, then oxygen has to settle out to pro- 
duce a nearly pure hyd rogen upper atmosphere (e.g. 
IKasting fc Pollackl 119831) . If oxygen is present in the 
evaporative flow at a high enough concentration, it will 
produce strong cooling and increase the opacity to the 
X-rays. This will severely suppress the gas tempera- 
ture, leading to a much lower evaporation rate, simi- 
lar to what occurs in a high metalicity protoplanetary 
discs (jErcolano fc C larke 20T(|), but more extreme. Now 
suppose all these conditions are satisfied and all hydro- 
gen in the water atmosphere is lost, this removes 1/8 of 
the atmosphere mass. However, since the original water 
atmosphere has a low scale height, such a removal can 
hardly change either the bulk size or the bulk density 
of the planet. One would therefore not be able to ex- 
plain the upper envelope in the observed planet radius, 
or the bimodal planet size distribution, or the correlation 
between planet density and X-ray exposure in terms of 
planetary evaporation. 

In contrast to atmosphere composition, we can be less 
certain about the core composition. For planets with 
hydrogen envelopes that are more than a percent in 
mass, the planet sizes are not sensitive to the core sizes 
(or equivalently, core composition). Planets that have 
been evaporated down to bare cores may be able to in- 
form us on the core composition, if the core masses are 
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Fig. 11. — The evaporation 'efficiency' (77) as a function of time. 
Here, we trace the evolution of four 20 Mjj planets, with core 
masses of 7.5 (solid), 10 (dashed), 12.5 (dotted) & 15 (dot-dashed) 
Mj, at a separation of ~0.03 AU from a sun-like star. All models 
lose their entire envelopes towards the end. 

known. A number of these objects show densities that are 
compatible with r o cky or iron/rock co mpositions (e.g., 
lHatzes et all 120111: iBatalha et all 120111 ) . More detailed 
investigations are necessary to ascertain the core compo- 
sitions. 
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Fig. 12. — The mass and radius evolution of a 20 Mgj planet with 
a 12.5 Mq core at a separation of ~ 0.05 AU undergoing standard 
Owen & Jackson (2012) evaporation (solid) -same as shown in Fig- 
ure [2} and pure EUV evaporation (dashed). The thin solid line 
represent the time at which the X-rays/EUV begin to decline. 



5.2. Improving the evaporation model 

As we have discussed previously, evaporation is key to 
the evolution of close in planets. Thus it is important to 
assess the role the assumed evaporation model plays in 
our conclusions. Most previous attempts to model the 
evolution of evaporating planets use a constant evapora- 
tion efficiency (77) , whic h is typically taken to be ~ 10% 
(e.g., lLopez et al.ll2012l ). We have argued that this effi- 
ciency depends on planet mass and radius, as well as on 
the X-ray flux. We further demonstrate this point here 
by showing how the efficiency changes as a planet evolves 
in Figure [TT1 Following four planets with the same initial 
mass of 20M® (but different core masses), we find that rj 
can decrease by a factor of 4 as the planets evolve from 
the early puffy stage to the later denser stage, though the 
variations are not strictly monotonic in time. The over- 
all decrease can be understood: as the planet evolves 
due to mass-loss and thermal contraction, the planet's 
density and the surface escape velocities increases with 
time. Therefore, it takes longer to accelerate the flow to 
the escape velocity. This leads to enhanced cooling and 
lower efficiencies. 

While the fixed ef fi ciency of ~ 1 0% a dopted by, e.g., 
Uackson et al.l (|2012D ; lLopez et al.l (|2012D does represent 
the approximate median value during the evolution, it 
can result in order unity inaccuracies in the integrated 
mass-loss. So any inferences about the initial planet 
structure should be taken with caution. 

5.2.1. Accuracy of the X-ray model 

The lOwen fe Jackson! (|2012D model for X -ray evapora- 
tion contain a number of assumptions that may impact 
our conclusions here. First, since only the soft < 1 — 2 
keV photons are responsible for heating, while the X-ray 
flux refers to the entire observed X-ray spectrum (0.1 — 10 
KeV). If the adopted X-ray spectrum - while based on 
observed spectra (see lOwen et alj l2010f )- is overly soft 
or overly hard, then the X-ray heating will be over- 



or under-estimated, respectively. Second, these calcu- 
lations assumed that, in the X-ray region, molecules are 
photo-dissociated or thermally dissociated by UV pho- 
tons. Therefore, Owen & Jackson neglected cooling as- 
sociated with molecular species. This is an important 
assumption which has not been vigorously tested. If in- 
stead the molecular species provide significant cooling 
(in the temperature range of 2000 — 5000K), the X-ray 
driven flow will not reach T > 5000K and may remain 
sub-sonic all the way out to the EUV ionization- front. 
This will markedly reduce the mass-loss rate. More de- 
tailed calculations are necessary to address this issue in 
the future. 

5.2.2. The role of EUV Evaporation 

We find that EUV evaporation contributes < 10% of 
the total mass-loss during the planet's evolution. The 
EUV portion of the stellar flux, though comparable in 
energy to that of the X-rays, are much less efficient at 
driving a high integrated mass-loss. This is because in the 
EUV flow region, the recombination of hydrogen and the 
subsequent cooling bleeds much of the energy to space. 
The radiation hydrodynamics is similar to that of a HII 
region where ionization is balanced by recombination and 
the tempera ture is controlled by the cooling thermostat 
to - 10 4 K ([Murray-Clay et al.ll2009h . 

We re-calculate the evolution of our 'standard' model 
(see Section l3~3|) undergoing EUV evaporation only, the 
results of which are shown in Figure [12] At high EUV 

fluxes (> 10 4 erg s" 1 ), the mass-loss scales as L ^ y / a 
(|Murrav-Clav et al.|[200llOwen fe Jackson! 120 12[ ) and is 
significantly lower than the values that apply to X-ray 
flow. This suppresses the mass-loss at early times by a 
factor of ~ 10, relative to the X-ray model. The inte- 
grated mass- loss is ~ 1 M®, as opposed to ~ 7 M® in 
the X-ray model. As such, EUV evaporation alone is un- 
able to sculpt the close-in planet population and cannot 
explain the observed features discussed in Section |4j 
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6. CONCLUSIONS 

In this work we have coupled the hydro dynamic evap- 
oration models of lOwen fe Jacksonl ()2012[ ) to the stellar 
evolution code mesa, in order to follow the mass and ra- 
dius evolution of low mass planets orbiting close to their 
stars. Evaporation, while having little effect on massive 
planets, can remove the entire Hydrogen envelopes for 
the hottest low-mass planets. In all cases we find that 
X-rays is the dominant sculpting force, and that most of 
the mass-loss occurs in the first few 100 Myrs when the 
stars were more chromospherically active and while the 
planets were still contracting thermally. 

Our main conclusions are as follows: 

1. Evaporation produces an upper envelope in planet 
radius as a function of separation. The location 
of this upper envelope depends on planet masses 
and the X-ray luminosity of the host stars. In par- 
ticular, M-dwarfs, having proportionally larger X- 
ray fluxes, should have stronger evaporative power 
than indicated by their low bolomctric luminosities. 

2. To closely reproduce the observed upper envelope 
in Kepler candidates, both around sun-like stars 
and around stars of other spectral types, we require 
that the most massive hot Neptunes should have 
total masses not much exceeding ~ 20 M@, and 
core masses roughly half of that. 

3. Very close-in planets can be stripped of their entire 
atmosphere. The boundary between complete loss 
and planets that can retain at least ~ 1% of their 
atmosphere lies at ~ 0.1 AU if the core masses are 
~ 6M@. At this distance (where most of the Ke- 
pler planets lie), a thinner envelope can not survive. 
So we expect a gap in planet size distribution - 
and this is observed in the Kepler catalogue, where 
there is a deficit of planets at ~ 2i? ffi straddling 
planets with high X-ray exposures and those with 
low X-ray exposures. Comparison with our mod- 
els suggests that most of the Kepler planets should 
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have core masses ~ 6Mq and should have primor- 
dial Hydrogen/Helium envelopes at least a percent 
in mass. Moreover, about half of the Kepler plan- 
ets belong to the high X-ray group and have been 
stripped down to naked cores. 

4. Evaporation naturally explains the observed corre- 
lation between planet density and separation. At 
closer separations, planets in general have higher 
densities. 

5. Combining all evidences that support evaporaton, 
we argue that Kepler planets were born with hydro- 
gen ^enveJo_p_es_j_jiot volatile-rich atmospheres (also 
sec IWu fc Lith wick 20ll. 

Looking ahead, we expect that the approach wc 
adopt here, coupling thermal evolution and evaporation, 
will perhaps be the only hope we have for recovering 
the initial structure of low-mass planets. With better 
determined stellar radii and hence more accurate planet 
radii, we may be able to retrieve the initial distribution 
of planet total masses and core masses. With more 
measurements of planet densities, it may be possible 
to reconstruct the histor ies of indiv i dual planets (as 
is done for Kepler- 11 by lLopez et al.l f2012l ). A larger 
sample of planets with measured masses and radii will 
allow us to place constraints on the core compositions, 
as well as initial planet entropy. It is at this point 
that we will begin to learn valuable information about 
the planet population at birth and make inferences 
about the planet forma ti on process (e.gllda fe LiiJl2005b 
IMordasini et"ai1l2012ai rbl:lH ansen fc Murravll2012ll2013| ). 
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